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A  class  of  statistics  for  the  two-sample  survival 
analysis  problem  is  introduced.  The  SW-statistic  may  be  written 
as  the  integral  of  a  weighted  difference  in  the  estimated 
survival  functions,  the  integral  being  with  respect  to  Lebesgue 
measure  on  time.  Since  the  integral  is  with  respect  to  real-time 
the  statistics  are  not  generalized  rank-statistics.  However, 
with  an  appropriate  choice  of  weight  function  they  are  non- 
parametric  in  the  sense  that  weak  convergence  is  guaranteed  for 
any  underlying  configuration  of  survival  and  censoring 
distributions . 

Asymptotic  distribution  theory  under  the  null  is  derived. 
Consistency  under  a  fixed  alternative  is  shown.  Efficacy 
expressions  under  natural  sequences  of  local  alternatives  are 
given  and  an  expression  for  the  most  efficient  weight  function  is 
developed.  The  asymptotic  efficiencies  of  some  specific  SW- 
statistics  under  the  proportional  hazards  alternative  are 
examined . 
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INTRODUCTION 

Various  nonparametric  test  procedures  for  the  two-sample 
censored  data  survival  analysis  problem  have  been  proposed. 

These  include  the  weighted  logrank  tests  (Gill  1980) ,  their 
supremum  versions  (Fleming,  Harrington  and  O'Sullivan  1986),  a 
generalized  Smirnov  test  (Fleming  et  al.  1980)  as  well  as  the 
orthogonal-score-statistics  procedure  of  Breslow,  Edler  and 
Berger  (1984).  These  procedures  are  based  on  generalized  rank 
statistics  and  hence  they  are  certainly  robust.  Indeed,  because 
they  are  based  on  ranks  they  are  invariant  under  monotone 
increasing  transformations  of  the  data.  However,  such  invariance 
is  not  always  a  desirable  feature  (O'Sullivan  1986) .  A 
difference  which  is  large  may  become  negligible  under  a  monotone 
increasing  transformation  of  the  data  (Figure  1) .  Hence  rank- 
statistics  will  not  distinguish  the  two  situations  of  Figure  1  in 
their  powers  and  yet  intuitively  in  applications  one  would  expect 
a  reasonable  test  procedure  to  detect  the  larger  difference  with 
more  power  than  the  smaller  difference.  A  class  of  non-rank- 
based,  nonparametric  test  procedures  for  this  problem  is  proposed 
here.  The  test  statistics  reduce  to  the  scaled  difference  of  two 
generalized  L-statistics  in  uncensored  data.  (The  generalization 
allows  random  as  well  as  deterministic  weight  functions  in  the 
usual  definition  of  L-statistics) . 

The  basic  set-up  and  notation  along  with  the  SW- 
statistics  themselves  are  introduced  in  section  1.  Section  2  is 
concerned  with  asymptotic  distribution  theory  for  the  numerator 
and  consistency  results  for  the  denominator  of  the  SW-statistic . 


In  section  3  these  results  are  exploited  to  yield  null 
distribution  theory,  consistency  under  a  fixed  alternative  and  an 
expression  for  the  efficacy  under  a  sequence  of  local 
alternatives.  Finally,  since  these  statistics  are  indexed  by  a 
set  of  weight  functions,  the  choice  of  an  efficient  weight 
function  is  discussed  in  section  4. 

1.  DEFINITIONS  AND  ASSUMPTIONS 

Let  (  (X  i  j,U  i  j)  ;  j  =  l,...,nj;  i  =  1,2}  be  2n  independent 
positive  random  variables  where  n  =  n  l  +  n  2.  For  i  =  1  and  2, 
suppose  X  ij  and  U  have  (right  continuous)  survival  functions 
S  j(  •)  and  C  j(  •)  respectively.  As  is  usual,  X  j  j,  U  ij  and  X  = 
min(Xij(  U  i  j)  will  be  termed  the  survival  time,  censoring  time 
and  observation  time  random  variables  respectively  for  individual 
j  in  group  i.  We  wish  to  test  the  null  hypothesis 
H  0:  S  :(  •)  =  S  2(  •)  based  on  the  data  (X  i  j,  I[X  t  j  <  U  j  j]  ? 
j  =  l,...n  j:  i  =  1,2},  where  I [A]  denotes  the  indicator  function 
of  the  set  A.  This  is  the  simple  two-sample  censored  data 
survival  analysis  problem  under  random  censorship. 

We  will  assume  throughout  that  S  •)  i  =  1,2  are 
continuous,  (though  not  necessarily  absolutely  continuous) .  On 
the  other  hand  we  allow  the  C  ,(  *)  i  =  1,2  tc  be  arbitrary 
monotone  decreasing  functions  on  (0,  ®)  with  C  ^0)  =  1.  Thus 
improper  censoring  distributions  are  allowed. 

Let  S  i(  •)  and  C  j(  •)  denote  the  Kaplan-Meier  estimators  of 
the  survival  functions  S  j(  •)  and  C  j(  •)  respectively.  Define 
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T  c  =  sup{t:  ( (S  t(t)  v  S  2(t)  )  A  C  :(t)  A  C  2(t)  )  >  0}  where  xvy 
denotes  max(x,y)  and  x  y  denotes  min(x,y) .  Both  S  1  and  S  2  are 
'well-estimated'  on  [0,TC]  in  the  sense  that  for  i  =  1  and  2 
E(sVt)I[t  <TC])  =  E(Si(t)I[t  <TC]). 

Thus  [  0 ,  T  c]  is  a  natural  interval  on  which  to  compare  S  1  and  S  2. 

The  interval  [0,TC]  includes  the  interval  [0,T]  where 
T  =  T  j  A  T  2  and  T  i  =  sup{t:  (S*  j(t)  a  c  4 ( t )  )  >  0}  .  Only 
information  on  [0,T]  is  used  by  the  generalized  rank  procedures 
mentioned  above.  Since  the  generalized  rank  procedures  in 
essence  compare  hazard  functions,  and  hazard  functions  in  both 
groups  can  only  be  estimated  when  the  risk  sets  in  both  groups 
are  non-empty  (i.e.  on  [0,T]),  then  [0,T]  is  the  natural  interval 
to  use  for  these  procedures.  However,  the  interval  [0,T  c]  is 
more  appropriate  for  procedures  which  compare  survival  functions. 

A 

Let  W(  •)  be  a  random  weight  function,  measurable  with 

respect  to  <r(X{y  I  ( X  t  j  <  U  t  j)  ;  j  =  l,...,ni;  i  =  1,2,}  where 

/> 

<r(A)  is  the  <j-algebra  generated  by  A.  (W(  •)  estimates  an 
underlying  deterministic  lebesgue  measurable  weight  function  W(  •) 
in  a  sense  to  be  defined  shortly) .  We  define  a  member  of  the  SW- 
class  of  statistics  to  be 


We  will  propose  test  procedures  for  the  two-sample  censored  data 
survival  analysis  problem  which  are  based  on  these  statistics. 
Note  that  by  a  simple  integration  by  parts,  SW  can  be  written  as 
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|*2 

1*1 

n  :n  2 

t  * 

t  A 

(/  W(v)  dv)  dS  2(t)  - 

(/  W(v)  dv)  dS  j(t) 

n  i  +  n  2 

0 

0 

t 

0  1 

0 

in  uncensored  data,  which  is  the  difference  of  two 

t  A 

L-statistics  if  W  is  deterministic. 


2.  DISTRIBUTION  THEORY  FOR  SW 


Initially  we  will  consider  the  component  of  SW  from  the 
i  1  h  group 

T 

'Jnj/  '  W(u)  (SA  i(u)  -  S  j(u)  )du. 

0 


For  ease  of  notation  the  subscript  i  will  be  dropped  for  now.  We 
seek  distribution  theory  not  only  under  the  null  but  also  under  a 
fixed  alternative  and  under  a  sequence  of  local  alternatives. 

Thus  we  allow  the  configurations  (Sn,Cn,Wn)  to  vary  with  n  the 
sample  size. 

Let  rn  be  the  population  parameter  corresponding  to  T. 

Thus 


rn  =  sup{ t :  S  n(t)  a  C  n(t)  >  0} 
and  for  some  limiting  survival  functions  S(  •)  and  C(  •) 

r  =  sup{t:  S (t)  a  C ( t )  >  0). 

We  require  that 

(Al)  rn  <  r  and  lim  Tn  =  r. 

n  -»® 


Indeed  in  most  cases  of  practical  interest  rn  =  t  V  n. 


Convergence  of  the  configurations  is  also  required  in  the  sense 
that  for  the  survival  functions  S(  •)  and  C(  •)/  and  the  weight 


function  W(  •)  , 


(A2) 


aim  |  I  w  n(t)  -  W(t)  |  |  (0<T!  =  o 
n  -»  a> 

lim  S  n(t)  =  S(t)  and  lim  C  n(t)  =  C(t)  V  t  £  [0,  t] 
n  4  od  n  ->  ® 


Assumptions  which  are  to  remain  in  effect  unless  otherwise 
specified  will  be  denoted  by  (A#)  throughout  the  text,  where  #  is 
an  integer. 

The  following  theorem,  which  is  basic  in  the  derivation 
of  asymptotic  distribution  theory  for  SW,  is  a  generalization  of 
Theorem  2.1  of  Gill  (1983).  The  proof  follows  the  steps  of 
Gill's  proof  quite  closely  and  hence  will  not  be  presented  here. 
Instead  the  reader  is  referred  to  Appendix  2  of  O- Sullivan 
(1986),  wherein  a  detailed  proof  is  given. 


Theorem  2 . 1 

Let  h  n,  n  =  1,2,...  and  h  be  non-negative,  non¬ 
increasing,  continuous,  bounded  functions  on  [0,  rn]  n  =  1,2,... 
and  [0,  r]  respectively.  Assume  that 

lim  h  n(t)  =  h(t)  0  <  t  <  r, 

n  ->  co 


1 T 

(h(v))  2 

- d  S(v) 

JO  (S  (V)  )  2(C(v)  )  - 


and 


1  im 
n  ■*  a> 


n 


’Tn 

(h  "(v)  )  2 

lim-  - ds  n(v) 

n  ->  o>  Jo  (S  n(v)  )  2(C  n(v)  )  " 


<J 
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If  Z  "(t)  =  -  n 


S  n(t) 


and  Z  ”(t)=  B  - 


0  (S  (v)  )  2(C  ( v)  ) 


for  a  standard  Brownian  motion  B,  then  on  D[0, 


(h  n(t  A  T)  Z  n(t  A  T)  )  =>  h(t  A  r)  Z  ®(t  a  T)  . 


If  also 


lim  J0  |d(hn(v)  -  h(v))  |  =  0  V  t  <  r 
n  -»<d 


then  on  D[  0 ,  <n] 


h  n(v)  dZ  n(v)  =>  h(v)dZ°°(v) 


t  A  T 

h 

Jo 


(2.2) 


Z  n(v)  dh  n(v)  => 


t  a  r 

r 

i 

Jo 


Z  “(vJdhCv) 


where,  by  definition 


h  (v)  dZ  “( V)  =  h  ( s )  Z  °°( s )  -  Z  ®(v)  dh  (v) 

0  Jo 


Conclusion  (2.2)  is  most  useful  in  determining  asymptotic 
distribution  theory  for  SW. 


Corollary  2 . 3 


Suppose 


(2.4)  lim  Wn(u)Sn(u)du 

n  -» ®  Jo 


W(u)  S  (u)  du  <  co 


|w  n 
Ju 


(v)  |S  n(v)dv) 


(2.5)  lim  - 


0  (S  n(u)  )  2(C  n(u)  ) 


dS  n(u) 


t  (Or 

|w 


(v)  |S(v)dv) 


0  (S  (u)  )  Z(C (u)  ) 


dS(u)  <  <D 


then  on  C[0,  ®] 


(2.6)  U  n(  •)  =  -4  n  W  n(u)  (S  (u)  -  S  "(u)  )du  =>  UOT(  ■) 

Jo 


where  U  ®(  •)  is  a  mean  zero  Gaussian  process  with  variance 


function 


<r2(t)  =  - 


'  t  A  T  ^2 

t  at  t  W  (v)  S  (v)  dv) 

,  u  , 

0  (S  (u)  )  2(C(U)  - 


dS  (u) 


Proof. 


Prn  Pt 

Let  h  n(u)  =  I  Wn(v)Sn(v)dv  and  h(u)  =  W(v)S(v)dv. 

t  A  T  Ju  Ju 

Then  U  n(t)  =  Zn(u)dhn(u).  If  the  weight  function  W  n(  •)  is 

Jo 

positive  then  the  conclusion  follows  from  (2.2)  of  Theorem  2.1 
ft  Ar 

with  U®(t)  =  Z  ®( u )  W  ( u )  S  ( u )  du .  The  result  also  holds,  if 

Jo 

W  n(  <)  is  not  necessarily  positive  as  is  shown  in  Theorem  2.3.1 
of  O'Sullivan  (1986)  . 

We  now  introduce  the  time-point  rc.  If  Ti  5  r,. 
i,i'  =  1  or  2 ,  i  t  i',  define 

f  r  j  i  if  S  j(  t  j)  =  0  and  C  4(  r ,)  >0 

Tc  = 

\  t  :  otherwise 


l v.' --v  > - v-v- v- v-v-v 


In  almost  all  cases  by  Lemma  A.l  tc  is  the  limiting  form  of  Tc 

P  P 

in  that  T  c  -»  rc.  Indeed  in  general  T  *  Tc  '  r  c  -*  t  •  r 

Since  U  n(  •)  is  tight  in  C[0,  ®] ,  Theorem  4.2  of  Billingsley 

(1968)  allows  replacement  of  (2.6)  with 


\  n 


T'T 


W  n(u)  (S  (u)  -  S  (u)  )  du  =>  U  «( 


Tc) 


For  the  remainder  of  this  paper,  unless  otherwise  specified  we 
assume  that 

(A3)  r  <  ®  and  I  W(  •)  io,t)  <  ® 

where  |  |  |  |  A  is  the  supremum  norm  over  the  set  A.  Certainly  the 

condition  that  r  <  ®  is  always  true  in  practice.  The  techniques 
used  in  this  paper  are  heavily  dependent  on  the  boundedness  of 
W (  •)  and  further  work  will  be  necessary  to  generalize  the  results 
to  include  random  weight  functions  in  SW  which  are  unbounded  in 
the  limit.  Note  however  that  Corollary  2.3  which  deals  with 
deterministic  weight  functions  is  general  enough  to  allow 
unbounded  deterministic  weight  functions. 

Now  we  replace  the  deterministic  weight  function  w  n(  •) 

A 

with  its  estimator  W(  •)  a  random  weight  function.  Suppose 


(A4) 


|W(  •)  -  W(  •) 


t  0  ,  T 


Tc> 


P 


n  -»  ® 


Lemma  2 . 7 


If 

(2.8) 


1 

j 

W(U)  -  w  n(u) 

0 


<  (Cn(u)) 


l  0  ,  T 


A 


T 


) 


(W(u)-W  n(u)  )  (S(u)-S  n(u)  )du  -»  0 


It  is  enough  to  show  that 


\  n  %  (C  n(u)  )  '  (S  (u)  —  S  n(u)  )  i  t  o  /  T  I  is  bounded  in  probability 
(henceforth  denoted  by  BP).  Let  h  n(v)  =  S  n(v)  -J  (C  n(v)  )  "  and  h(v) 


=  S  (v)  ^C(v)  )  ‘  .  Then 


(h(v)j 


-1  im 
n  -*  ®  Jo 


h  n(v) 


dS(v)  =  1  -  S  (  t)  <  oo,  and 


dS  n(v)  =  1  im  (1  -  S  n(  r  °)  )  =  1  -  S(r). 


Thus  the  conditions  of  Theorem  2.1  which  guarantee  convergence  of 
h  n(  •  T)Zn(  •  '  T)  are  satisfied,  except  that  h  and  hn  are  not 
necessarily  continuous  but  rather  are  left  continuous.  Lemma 
2. A. 6  of  O’Sullivan  (1986)  shows  that  under  these  conditions 
h  n(  .  '  T)  Z  n(  •  A  T)  is  stable  (or  ’tight’)  at  the  end-point  r, 
even  if  h  and  h  n  are  only  left  continuous.  That  is  to  say,  given 
5  >  o  and  t  >  0  3  t*  <  r  and  n*  such  that  V  n  in* 

P[  I  I  *  ^Cn(  •))  '  (S(  •)  -  S  n(  •))  I  I  tt*,T]  >  f]  <  «• 

. —  A 

!  'h  (S  (u)  -  S  n(u)  )  |  I  t  o ,  t  *  ]  BP  as  a  consequence  of  Theorem  2.1. 

Thus  !  -J[C  n(  •) )  ’(S(  •)  -  S  n(  •) )  |  |  [0,T]  is  B.P.  as  required. 


The  conditions  on  W(  •)  implied  by  (2.4),  (2.5)  and  (2.8) 
certainly  guarantee  the  desired  weak  convergence.  However  more 
direct  and  interpretable  conditions  on  W(  •)  are  desirable  to 
facilitate  the  generation  of  SW  statistics  which  are  stable 

A 

asymptotically.  We  seek  conditions  on  W(  •)  which  will  yield  SW 


r*.  *-  *►*»>*./. 

J  J  J  ^ 
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asymptotically  stable  for  all  choices  of  underlying  survival  and 
censoring  configurations  (S  n(  •)  ,  C  n(  •) )  .  In  this  sense  the 
proposed  tests  will  be  non-parametric. 

Definition 

For  two  sequences  of  stochastic  processes  X  n(  •)  and  Y  n(  •) 
defined  on  [0,  ®)  ,  by  X  n(  •)  =  0(Y  n(  •) )  we  mean  that  3  a  constant  V 
such  that,  given  «  >  0  3  an  integer  n  e  with 

P[  |X  n(t)  |  <  T|Y  n(t)  |  V  t  >0,  IfnJiJ  >1-  f 

Note  that  if  X  n(  •)  and  Y  n(  •)  are  deterministic  functions,  then 
the  above  definition  reduces  to  the  existence  of  a  constant  T 
such  that  |Xn(t)  |  <  r|Yn(t)  |  Vt  >0,  V  n. 

Theorem  2 . 9 
Suppose 

W  n(  •)  =  0(  (  (C  "(  •  *  t")  )  ~)  1  '  2)  ,  then  on  C[0,  «d] 
ft  A  T 

in  J  W  n(u)  (S  (u)  -  S  n(u)  )  du  =>  U®(t) 

Under  the  stronger  condition  that 


1  *  **«•  *  -  *  *,**4  •*«_  *’«  f ^  •**  * «*,  «r  «**  ■'»  * 


( A5 )  3  3  >  0  with 

W  n(  •)  =  0(((C  n(  •  ‘  r  "))“)■  5  **) 


and 


then 


W(  •}  =  0((C  '(•))  ' 5f5)  , 


IT  '  T  c  *  r  c 


n  n 


(W(u)  -  Wn(u))(S(u)  -  S  n(  u )  )  du  ^  0. 


JO 


v  A 


’•  '»  *»  %  “»  %  *■ .  %  %  jf 


•»  v ’.■v/  ‘A*v.«’vv/ 

4.  <4.4  <4.  •-  ■O  ' 


Proof 


By  (A3),  (2.4)  is  satisfied.  For  some  constant  f, 


»  T 

(  r 

'  7 

7 

I  |W(v)  jS  ( v)  dv 

2 

/  S (v) dv 

u 

ds (u)  <  -r2 

t 

u 

0 

(S  (u)  )  2(C  (u) ) 

0 

(S (v) )  2 

dS  (v) 


=  r2  ff2(S(  *)I[0  <  *  <  T])  <  m 

where  <r2(S(  *)I[0  -  *  <  t]  )  denotes  the  variance  calculated  from 
the  survival  function  S  (  •)  I  [  0  <  •  <  r]  and  the  equality  follows 
from  Lemma  A. 2  of  the  appendix.  Define 

'  7  ) 

!  |W(v)  |S(v)dv  2 

^  u 

g(u)  =  - I[0  <  u  <  r] 

(S(u)  2(C(u) )  “ 


'  j  n  ' 

;  |w  n( v)  |s  n(v)  dv  2 

.  U  / 

g  n(u)  ■  - 1(0  <u  <  Tn]. 

(S  n(u))  2(C  n(u) )  ' 

Also  let  Yn=XnI[0  <Xn<  rn]  and  Y  =  XI  [  0  <  X  <  t]  where  X  n 

and  X  have  survival  functions  S  n(  •)  and  S(  •)  respectively.  The 
survival  functions  for  Yn  and  Y  are  (Sn(t)-Sn(  r"))  and 

d 

(S  ( t)  -S  (  r) )  for  t  >  0  respectively.  Since  Y  n  -»  Y,  g  is 
continuous  at  0,  g  has  at  most  a  countable  number  of  disconti¬ 
nuities  on  (0,  ®)  and  S(  •)  is  continuous  on  (0,  as),  then 
d 

g(Y  n)  ->  g(Y)  .  Given  e  >0  3  n  0  and  t  0  <  t  such  that  V  n  >  n  0, 

P[Y  n  <  t  0]  >  1-  f.  Also,  lim  |  |g  n(x)  -  g  (x)  I  I  t  0  ,  t  ]  =  0.  Hence 
p  n^®  d 

9  n(Y  n)  ~9(Y  n)  0  and  so  9  "(Y  n)  +  9(Y )  •  Boundedness  of  g  n(Y  n) 

uniformly  in  n  yields  lim  E(g  n(Y„))=E(g(Y))  and  hence  (2.5). 

n  ->  ® 

Finally  we  need  to  verify  that  the  stronger  condition 
(A5)  yields  (2.3).  (2.8)  is  a  trivial  consequence  of  (A4)  unless 


a 


C  "(  r  A  tc)  =0,  in  which  case  we  note  that  r  =  rc.  Thus  suppose 
C  "(  r)  =0  and  r  =  re.  Given  e  >0  3  t  0<  r  such  that  0  <  (C  "(t  q)  )  a 

<  t/2. 


W(u)  -  W  n(u) 

W(U)  -  W  n(u) 

< 

4  (C  n(u)  )  - 

f  0  ,  T  A  T  e  ]  4c  n(U)  ) 

t0,t0  A  TaT{] 

W(u)  -  W  n(u) 

4*  - 

4C  "(U)  )  t  t  Q  ,  T  A  T  £] 

The  first  term  on  the  right  hand  side  converges  to  0  in 
probability  as  a  consequence  of  (A4).  For  the  second  term  3 
constants  Tj  and  f2  such  that 


W(u)  -  W  n(u) 

4c  n(u) )  - 


t  t  o  ,  T  A  T  e] 


C“(u) 

<  r2  -  (c-(to))!+  r:((cn(t  o) )  -)  1 

(C  n(u)  )  -  {t  /TJ 

<  r2  +  r i  f 

with  probability  at  least  1  -  3  €  if  n  is  large  enough,  using 
Theorem  3.2.1  of  Gill  (1980)  for  the  first  component.  Thus 
I  |W(u)  -  W  n(u)  I  I  p 


4c  n(u)) 


l  o  . 


Ct0,T  a  T  , 


The  final  step  in  the  derivation  of  asymptotic  distribution 
theory  is  to  show  that  the  contribution  to  SW  over 
(T  ‘  T  c  ■'  rc,  T  c)  is  asymptotically  negligible. 

i  t  P  r  n  a  T 

C  W(u)  (S(u) -S  "(u)  )du  =  +  '  W(u)  (S(u) -S  n(u)  )du 

JT'Tc'tc  J  T  *Tt  '  rt 


V,'-.-',  >Vvv'-V-'. .A  .  V-V-Y- 


T 


c 


+ 


W (u)  (S  (u)  -S  n(u)  )du. 

7n  A  T  . 
c 


We  will  treat  these  two  components  separately. 


Lemma  2 . 10 
If 

( A6)  -^i(T  -  Tn)  is  BP,  then 

i  r n^I>  e 

-f  W(u)  (S(u)  -  S  n(u)  )  du  lo 

J  T^T  c  A  T  c 

Proof. 

For  some  constant  T  with  high  probability 
i  7n*r  c 

Ap  W(u)  (S(u)  -  S  n(u)  )du 

J  T  A  T  c  A  7  e 

<  r  :  'S(u)  -  S  "(u)  |  |  ( 0  ,  r ,  |  aF(  Tn  -  T)  I  and 
S  (u)  -  Sn(u)  !!io,r)  0  by  Gill  (1980)  Theorem  4.1.1.  [] 

A 

Remark.  The  uniform  consistency  of  S  (  •)  on  [0,  r)  is  stated  but 
not  proved  in  the  reference  cited  above.  A  slight  modification 
of  the  proof  of  Theorem  1,  page  304,  Shorack  and  Wellner  (1986) 
yields  the  result.  An  alternative  proof  may  be  obtained  from  the 
author  and  will  be  included  in  a  forthcoming  monograph  by 
Fleming,  T.R.  and  Harringdon,  D.P.  Condition  (A6)  holds  very 
generally  at  least  under  a  fixed  conf iguat ion ,  as  is  shown  in 
O'Sullivan  (1986)  using  some  results  from  extreme-value 

Lemma  2^11 

( A7 )  If  7t  =  r  2  assume  that  -*i(  r"  -  r^)  is  bounded,  then 


W(u)(S(u)  -  S  n(u) ) du  -»  0 

J  rn  at  , 


Proof , 


_  » T  c  _n 

(2.12)  -ft  W(u)  (S(u)  -  S  n(u)  )  du  =  -  -fr 

r  n  A  T  c  t 


C  A  n 

du  =  - -fr  W(u)S  (u)du. 


Tn  AT  c 


Under  a  fixed  configuration,  either  S(  t)=0  so  that  (2.12)  is 
identically  zero,  or  S(  r)  > 0.  In  this  case 

lim  P[SA(T)  =  0]  =  0,  lim  P[C(T)  =  0]  =  1  and  lim  P[T  c  >  r]  =  0 
n-*<»  n  -»<x>  n od 

so  the  result  holds. 

Under  a  sequence  of  local  configurations,  if  Tj  <  r3_j 

,  P 

for  l  =  1  or  2,  then  S  (  r  4)  >0  and  Tc  -»  tc  =  Tj  by  Lemma  A.l. 

Hence  lim  P[T  c  =  T  j  <  rn  a  rn]  =  1  and  the  result  holds.  If 
n  -»a>  12 

t  !  =  t  2  ~  T  it  is  enough  to  consider  rn  <  t  V  n. 

IT  c 

-\J  n  W(u)  S  n(u)  du  <  ^i(  rn  -  r ")  S  n(  t")  |  |w  (u)  |  | 

JrnATc  1  2  (rn,Tc) 

If  S(  t)  =0  then  the  result  is  proven.  Otherwise  S(  t)  >0  and 

C  n(  t°)  =0  for  n  large  enough.  In  this  case  C(  rn)  I  0  since 


|C(u)  -  C  n(u)  |  |  t0#T)  I  o,  and  |  |W(u) 


[  rn,T 


We  are  finally  in  a  position  to  give  asymptotic  distri¬ 
bution  theory  for  SW.  The  natural  condition 


(A8)  lim  - =  P{  >  0 

n  -» oo  n 

is  imposed.  Note,  that  the  assumptions  listed  so  far  were  given 
without  regard  to  a  group  subscript  i.  We  assume  that  the  condi 
tions  hold  for  both  groups  i  =  1,2. 


y 


Lemma  2 . 13 


n  in  2 


i  +  n  2 


C  A  A  A 


(W(S  ,  -  S  2)  -  W(S  n  -  S  n)  )du  ->  N(0,  <r‘  ) 


l  2 


where 


Proof . 


-  V 


W(v) S (v) dv) 


0  (S  [(u)  )  2(C  i(u)  ) 


dS  i(u) 


By  Theorem  2 . 9  and  independence 


n  ln  2  f  ft  AT  x  A  It  AT  2 

-  W  n(u)  (S  t(u)-S  "(u)  )du  -  W  "(u)  (S  2(u)-s  "(u)  ) 

n  i  +  n  2  IJo  1  Jo  2 


=>  ^  p 2  U  j(t)  +  ^  Pi  u  2(t)  on  C [ 0 ,  €»]  . 


That  T  c  A  t c  ■*  rc,  the  sequence  is  tight  in  C[0,  o>]  and  Theorem 


4.2  of  Billingsley  1968  together  imply 


n  p  2 


l  c  '  e 


T  ^T  c a  T  c 


•  -  Wn(u)  (S  i(u)-S  "(u))du  -  W  n(u)  (S  2(u)  -S  "(u)  )du 

n  j  +  n  2  l«0  1  Jo  2  ) 

-*  ■'J  P2  U  i  (  T  c)  +  -J  Pi  U2(Tc) 

=  P2  ^  l(  T  l)  +  a]  Pi  U  2(  T2)  • 

Theorem  2.9,  Lemma  2.10  and  Lemma  2.11  then  yield  the  required 
result. 


Having  established  the  asymptotic  behaviour  of  the 
numerator  the  next  task  is  to  propose  consistent  estimators  of 


</•  f  v  A  .  v.v.v.  -- v.\  . 


'A  ' -• . v  . 


•V'V’.'.V'AV.V-V. 


n  i  -  1 


luu  w ( v)  S  i(v)dv)  1 


S  }(u)  S  i(u)C  |  ( u) 


IS  j(u)  . 


A  A  C 

where  pi  =  n  /n.  <rup  is  termed  the  unpooled  variance  estimator. 

A  P*  C  A  A  — 

In  uncensored  data  with  W(  •)  =1,  j  W(u)S  4(u)du  =  X  if  the  sample 
mean  from  the  ith  group  and  Lemma  A. 2  yields 


n  i  ~  1 


i  {  i  A  A 

Uu  W  (v)  S  i(v) 


0  S  i(u)S  i(u)C  i(u) 


dS  i(v) 


^  (x  i  j  -  xi.) 
£—>  n  ;  -  1 


the  usual  sample  variance  of  survival  times  in  the  ith  group. 

A  2 

Thus  <rup  is  the  natural  estimator  of  the  variance  of  SW  at  least 
in  this  classical  case. 

Under  the  null  hypothesis  it  is  tempting  to  substitute 
S pf  the  Kaplan-Meier  estimator  calculated  using  the  two  samples 

A  A  2 

pooled,  for  S  j  in  <rup  .  The  resulting 


,2  _  _ 


n  -  1 


(£v 


(v)S  p( v) dv 


dS  D(u) 


i=l  JO  S  p(u)S  p(u)C  j(u) 

is  termed  the  pooled  variance  estimator.  Indeed  simulation 

A  2 

results  suggest  that  <rp  is  a  better  standardization  factor  for  SW 
/  2 

than  <rup  under  the  null  in  many  cases  in  small  sample.  We  defer 
discussion  of  small  sample  properties  to  another  paper  however. 


Theorem  2 . 14 


a  2  p  2 

p  J)V 


/  v  .*  ‘y  v  .*  v  ',■«  • v  .*  •  *  ^  •  v,  -  ,  - .  - « * .  • .  -  ,  ■  .  -  »*- .  %  * ,% 

-  A  /-  .'-.'■.'•'.'•'a''" .N.'IV 


!  8 


Under  H  0:  S  t  = 


S  2  =  S  V  n 


A  2  p  2 


The  proof  of  Theorem  2.14  is  somewhat  technical  and  is  given  in 


O'Sullivan  (1986).  Consistency  of  <rp  is  proven  only  under  the 


null  hypothesis  since  general  results  on  the  consistency  of  the 


Kaplan-Meier  estimator  are  known  only  under  the  random  censorship 


model.  Unless  the  survival  distributions  in  the  two  groups  are 


equal  or  the  censoring  distributions  in  the  two  groups  are  equal 


then  the  random  censorship  model  does  not  hold  in  the  pooled 


sample.  Consistency  of  the  Kaplan-Meier  is  basic  in  deriving 


consistency  of  the  variance  estimators.  Thus,  consistency  of  <rp 


is  not  proven  in  general  but  is  shown  only  under  the  null. 


Although  consistency  can  also  be  shown  for  arbitrary  config¬ 


urations  under  the  assumption  of  equal  censoring  distributions 


n  n 

(C  j  =  C  2  V  n) ,  this  assumption  is  very  restrictive  and  we  do  not 


consider  it  to  be  an  important  case. 


3.  NULL  THEORY,  CONSISTENCY  AND  EFFICACIES 


The  purpose  of  this  section  is  to  summarize  the  results 


of  Section  3  as  they  pertain  to  null  distribution  theory, 


consistency  and  asymptotic  behaviour  under  a  sequence  of  local 


alternatives.  The  following  assumptions  are  made,  for 


i  =  1  and  2 


<  r  j  V  n,  lim  rn  = 


>  A  .  •  -  •  .  *  .  ■  .  ■  A  .  - 


V.  *\  -  .  s.  V,  . 

,  ',r.  *"r  . 


I -  n  ,  ,  . 

(A6)  -fi  (T  j  -  Tj)  is  bounded  in  probability  , 

—  n  n 

( A7 )  If  Tx  =  r2  then  'O  1  (  T1  ~  T  2)  is  bounded 

n  i 

( A8 )  lim  - =  pj  >  0. 

n  ®  n  j  +  n  2 

Under  a  fixed  configuration  (Al) ,  (A2)  and  (A7)  are 
redundant,  and  (A6)  holds  except  in  pathological  cases.  In 
practice  Tj  <  ®  and  (A5)  implies  that  W(  •)  is  bounded.  (A8)  also 
holds  in  general  since  the  design  of  the  experiment  will  usually 
ensure  that  (n  /n)  is  a  fixed  positive  fraction.  Thus,  only  (A4) 
and  (A5)  need  be  of  real  concern  for  the  purposes  of  null 
distribution  theory  and  for  consistency. 


For  a  sequence  of  local  alternatives  it  is  most  natural 


to  consider  the  censoring  distributions  as  fixed  and  the 

alternatives  are  specified  in  terms  of  the  convergent  sequences 

of  survival  functions.  In  most  cases  of  practical  interest  for 

i=l  and  2,  C  ^  r  j)  =  0  and  S  j(  r  >  0,  i.e.  there  is  a  positive 

probability  of  surviving  past  the  end  of  the  study.  Thus  even 

under  a  sequence  of  configurations  (Al)  and  (A7)  are  redundant  in 

n 

most  cases  of  interest  since  Tj  =  t,  V  n,  i=l  and  2.  As 
discussed  above  for  a  fixed  configuration  (A3)  and  (A8)  are 
generally  satisfied  in  practice.  Thus  for  a  sequence  of  local 
alternatives  only  (A2),  (A4),  (A5)  and  (A6)  are  of  any  real 

concern  in  general. 

Let  SW  up  =  SW/  %'<r2  and  SW  p  =  SW/  ^a2 

up  p 

Theorem  3 .  1 

(i)  Under  the  null  hypothesis 

d  d 

SW  up  -»  N(0, 1)  ,  and  SW  p  ■»  N(0, 1)  . 

(ii)  A  two-sided  test  based  on  SW  u  p  or  SW p  is  consistent  against 

I  PTC  ! 

any  fixed  alternative  such  that  |  J  0  W (u)  (S  :(u)  -  s2(u))du|  >  0 

(iii)  If  under  a  sequence  of  local  alternatives,  for  some  bounded 
function  D(  •)  on  [0,  r  c)  , 

i  n  ln  2  n  n 

-  (S  x(u)  -  S  2(u))  ->  D(u) 

n  j  +  n  2 

uniformly  on  [0,  r  ),  then 


'■■■II  l»u 


I 


SW  u  P  -»  N 


fc 


W (u) D(u) du 


,  1 


S  V 


Proof 

(i)  This  is  a  trivial  consequence  of  Lemma  2.13  and  Theorem  2.14 

a  2  p  2  ,2 

(li)  <rup  a  <rsw  <  co  so  that  ffup  is  BP.  Also, 

C  r  /  2  ,  , 

W(v)S  p(v)  dv  !  r  P  2  P  l 


P  T  c  /• 


2  _  _ 


V  dU 


0  Sp(u)Sp(u) 
F  c 


C  x(u) 


C  j(u) 


d  S  _(u) 


<  -  r 


\2 


\  uU 


S  p(v)dv 


JO  S  p(u)  S  p( u) 

T  with  probability  at  least  1-  £  V  n  <  n  f, 

S  (  Tj  v  t2)  2  <  cb. 


dS  p(u)  for  some  constant 


Thus  it  remains  to  show  that  |SW  |  -» 


IT 


C  A 


Pt, 


W(u)  (S  i ( u )  -  S  2(u)  )  du  - 


JO 


W(u)  (S  j(u)  -  S  2(u))du 


Fca 

(W(U)-W(U)  )  (S  X(u)-S  2(u)  )du 

+ 

W(u)  (S  i ( u )  -S  2(u)  )  du 

Jo 

Jt  e  atc 

Pt, 


W(u)  (S  j(u)  -  S  2(u)  )du 


|  JT  CA  Tc 

The  first  term  converges  in  probability  to  zero  by  (A3)  and  (A4). 

Lemma  A.l  and  boundedness  of  W(  •)  prove  that  the  last  term 

.  .  .  P 

converges  to  zero  in  probability.  Similarly  since  T  5  -»  t  c  unless 
T  i  <  T3  -  i  with  S  j  (  r  j)  =0  and  C  j(  r  j)  =0  (Lemma  A.l),  then  the 


Si 

t| 


A 

I 

1 


:A 


A 

9* 


w*, 

v.*, 

£ 

I 

rrt 


■i 


•  •  -SS'S  •  •  •  “  v"  v"  \  ;■  • ,  -y\ 

.  •  *  »  «  *  .  -  *  »  »  -  .  .  ■  »  ,  ««•  «  *  *  .1  >  "  •  •  -  *  .  "  -■  /•  ,  f  "  V* .  »  w  w  »  w  -«  <  »  -  _  •* „ 

A  a  f»  it  1 1  Wh  mrf  n  fti*4  i  At  i  ift ■  »•*  et*  tCm  ifli  tCm  if d  ffl&jwV L  *j(  jlA 


p[C  i(  rt)  =0]  •*  1. 


n  ln  2 


fT  cA  T  c 


n  1  +  n  2  »0 


(W(u)  *  W(u))(s"(u)  -  S  2(u)  )du  I  0 


r 


n  ln  2 


since 


|S  "(u)  -  S  2(u)  |  |  t  0  t  T  )  is  bounded  and  ( A4 ) 


S  n  j  +  n  2 

holds.  Thus  (4.2)  is  asymptotically  equivalent  to 


r“ - 

TcATc 

n  ln  2 

n  n 

W(u)  (S  t(u)  -S  2(U)  )  du  - 

N  n  j  +  n  j* 

0 

W  (u)  D  (u)  du . 


e  p 

=  -  W(u)  D  (u)  du  ■*  0 

Jt  e  '  r  c 

fPTc  )  2 

j  W(u) D (u) du  I  1 

\  Jo 

Note  that  e(SWup)  =  -  is  termed  the  efficacy  of 

,!1V 

SW  u  p  under  the  sequence  of  alternatives  specified  by  D(  •) . 

4.  EFFICIENT  WEIGHT  FUNCTIONS 

The  objective  in  this  section  is  to  find  a  weight 
function  Wopt(  •)  which  maximizes  e(SWup)  over  all  possible  weight 
functions  W(  •)  for  a  fixed  S(  ■)  ,  C  x(  •)  and  C  2(  •)  and  for  local 
alternatives  specified  by  D(  •)  .  Since  e(SWup)  is  invariant  under 
scalar  multiplication  of  W(  •)  the  constraint  is  imposed 


«  -.lIvSvS 


(4.1)  1  =  <r 


2  _  _ 


I  W 

du 


(v) S(v)dv  I  4 


0  S  (u)  S  “(u) 


P  2  Pi 


C  i(u)  C  2(u) 


dS (u) . 


Direct  maximization  of  (  /  W(u)D(u)du)  2  with  respect  to  W(  •) 

0 

subject  to  <r 2  =  1  is  a  difficult  problem  in  calculus  of  varia- 

$  w 

tions  and  indeed  the  Euler-Lagrange  equations  are  in  general  not 

solvable.  An  alternative  approach  is  to  transform  the  problem  to 

uo>  Tc 

the  maximization  of  e(SWup)  with  respect  to  h(  •)  =  f  _  W(u)S(u)du, 
subject  to  the  additional  constraint  that  h(  tc)  =  0.  Some 
regularity  will  be  required.  The  optimal  h(  •)  is  denoted  by 

h  ,pt(  ’)  • 


Lemma  4 . 2 

D(u) 

If  lim  -  <  ®  ,  S  (  •)  and  D(  ’)  are  differentiable  and 

U  *  7C  S (U) 


dS(u)  *  0  V  u  <[0,tc),  then 
du 


hopt(u)  = 


S  2(u)  C  !(u)C  2(u)  du  \S(u) 


,  U  <  T  c, 


P  ]C  t(u)  +  Pf  2(u)  d_  (S  (u)  ) 

du 


h  oPt(  Tc)  -  0. 


Proof 


Let  V (u)  =  D(u)/S(u).  Since  D(0)  =  0,  V(0)  =  0  and 


integration  by  parts  yields 


W(u) D(u)du  =  -  V(u) dh (u)  =  h(u)V'(u)du 

Jo  Jo  Jo 


(  pf  j(u)  +  pfi  2(u)  \ 

where  V'  (u)  =_d_  (V(u)  )  .  Let  A(u)=  -  d  (S(u)  )  - 

du  du  S  2(u)C  "(u)C '(u)  I 


Then 

f  Ptc  ]2  ( prc  V'  (u)  ]2 

I  W(u)D(u)du  =  -  ^(u)  h(u)du 

iJo  /  l  Jo  •A(u)  / 

(  prc  (V*  (u) )  2  'j  (  prc 
<  - du  A(u)h2(u)du 

Jo  A(u)  IJo 

Tc  (V* (U) )  2 

- du 

0  A  (u) 

where  the  first  inequality  is  an  application  of  the  Cauchy- 
S\-hwartz  inequality  and  the  second  equality  follows  from  the 
constraint  (4.1).  The  inequality  is  an  equality  if  and  only  if 
V  (u) 

-  a  \A(u)  h  (u)  ,  u  i  [0,  Tc). 

*(u) 


If  lim  hopt(u)  #0  then  the  solution  requires  h  0  p  t  to 
u  *  r  c 

have  a  mass  point  at  rc.  If  censoring  is  continuous  at  rc  or 

S(rc)  =  o  then  certainly  lim  hopt  (u)  =  0  if  V'  (u)  is  bounded 

u  *  t  e 

on  some  (  tc  -  t,  r  c)  ,  <  >  0.  However  if  hopt  is  not  continuous 

at  t  c  then  an  optimal  weight  function  W(  •)  will  not  exist.  Under 
the  assumption  that  hopt(  •)  is  differentiable  then  the  optimal 
weight  function  Wopt  which  maximizes  e(SW)  is  given  by 


2b 


1 


S  2(t)  C  j(t)  C  2(t) 


(D(t)/S(t  ) 


S(t)  dt 


|  d  S(t) 
Pf  t(t)  +  Pf  2(t)  1  dt 


Even  if  an  optimal  weight  function  exists,  Wopt(  •)  need 
not  satisfy  the  reguxarity  conditions  of  section  3.  For  example 
Wopt(  •)  need  be  bounded  on  [0,  tc).  However  it  is  clear  from 
Corollary  2.3  that  the  weaker  conditions  (2.4)  and  (2.5)  will 
yield  weak  convergence  for 


n  in  2 


n  +  n 


1AT  2 


W0pt(u)  ((S  t(u)  -S  2(u)  )  -  (S  i(u)  -S  j(u))du 


since  Wopt(  •)  is  deterministic  and  T  2  ->  i-  Tt  a  Tj.  Under  a 
sequence  of  local  alternatives  r=  tc,  so  that  it's  asymptotic 
efficacy  is  given  by  e(SWuP).  Thus  Wopt(  •)  I  [  •  <  T  r'T  2]  truly 
yields  an  optimal  SW-statistic  for  a  given  S,  Cx  and  C2  and 
sequence  of  local  alternatives  specified  by  D,  under  these  mild 
conditions. 

Such  optimal  weight  functions  do  not  however  yield  non- 
parametric  statistics,  in  the  sense  that  for  a  fixed  wopt,  weak 
convergence  of  SW  is  not  guaranteed  for  all  choices  of  (S,  C 
C  2) .  It  is  still  informative  to  calculate  the  optimal  weight 
function  and  indeed  the  behaviour  of  a  particular  non-parametric 
SW-statistic  can  be  explained  in  part  by  a  comparison  of  the 
nonparametric  weight  function  with  the  optimal  one. 

A  nonparametric  weight  function  which  we  propose  for 
general  use  is  given  by 


mm 


>v.\vNv.v,y  -V- . 


c  x(  -)c  2(  •)  j 

w  c(  •)  =  -  • 

A  A  ~  A  A  " 

Pf  !(  *)  +  P  2(  0 

In  uncensored  data  W  c(  •)  =  1  so  that  the  corresponding  SW 
procedure  is  essentially  a  generalization  of  the  z-test  to 
censored  data.  In  the  following  example  we  compare  the  behaviour 
of  this  SW  statistic,  and  the  optimal  SW  statistic  to  the  locally 
most  powerful  test  under  the  proportional  hazards  alternative 
with  equal  censoring  distributions,  namely  the  Logrank  test. 


Example  4 . 3 

Let  Sa(t)  =  (S0(t))a+1,  a.  >  -1.  A  sequence  of  local  pro¬ 
portional  hazards  alternatives  is  given  by  S  "(  *)  =  S  0(  ■)  , 

S  ”(  •)  =  S  an(  •)  with  a n  =  a  n  x  +  n  2  for  some  constant  a-  For 

■m  n  5n  2 

any  t  such  that  S  0(t)  >  0  then 

n  tn  2  n  .  ♦  n  2 

-  (S  n(t)  -  S  n(t)  )  =  -  a  S  0(t)  log  S  0(t)  *  0  - - 

^  j  ♦  n  j  1  2  ,  n  ;n  2 

-  -  a  S  0(t)  logs  0(t)  =  D(t)  .  If  S(t)  logs  (t)  ,  <  ®  then  the 

convergence  is  uniform  on  the  set  A. 

If  the  underlying  survival  distributions  are  Weibull, 

.  k  k  .  .  k 

S  ,(t)  =  e  ,  k  >  0,  then  D(t)  =  a  t  ‘e  and  the  convergence  is 
uniform  on  [  0 ,  u)  .  Simple  algebra  yields  Wop?(t)  =  t  1 "  *  in 
uncensored  data  and  Wopt(t)  *  1  +  k  ( t  ,  -  t)tk_1,  t  *  ;  0 ,  t  y  , 
under  equal  uniform  censoring  on  [0,t  '  The  asymptotic  efficacy 


of  the  (most  efficient)  Logrank  statistic  is  given  by 


e  ( Lgk) 


pr  C  ! ( u )  C  j ( u ) 

- - - - - dS  0(u) 

Jo  PlC  -(U)  +  f>£  *(u) 

1  2 


(Gill  1980) , 


so  that  the  asymptotic  efficiencies  (A.E.)  of  SW c  and  SW 0 p t  are 
given  by  e  ( SW  c) /e  (  L  g  k)  and  e  (SW  0  p  t) /e  ( L  g  k)  respectively.  These 
were  calculated  under  a  variety  of  configurations.  The  results 
are  tabulated  in  Table  1. 


Table  1 

Asymptotic  Relative  Efficiencies  under  Weibull  Proportional 
Hazards  Alternatives 
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AE  ( W  c) 
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1) 

87 
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3) 
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1) 
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(  3  , 

3) 
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(  3  , 

J 

.  84 

1 

.  00 

.  denotes  the  uncensored  case 


In  all  cases  examined  SW 0  f  is  fully  efficient  and  SW  . 
maintains  high  efficiency  across  this  broad  range  of  survival 


configurations.  In  the  classical  exponential  scale  family  with 
no  censoring  Wopt(  •)  =  1  *  K  i*)  and  the  z-test  is  in  fact 

optimal.  For  values  of  k  greater  than  unity  the  z-test  loses 
efficiency  in  uncensored  data.  The  optimal  weight  function 
emphasizes  later  differences  in  the  survival  functions  more  than 
early  differences  and  behaves  less  like  W  c(  •)  =  1  as  k  increases. 
In  censored  data  the  AE  of  SW  c  also  decreases  as  k  increases, 
again  because  Wc  mimics  the  behaviour  of  Wopt  more  closely  for 
small  values  of  k  than  it  does  for  large  values. 


CONCLUDING  REMARKS 

From  the  point  of  view  of  the  applied  statistician  the 
SW-stat ist ics  are  an  intuitively  appealing  class  of  statistics 
for  the  two-sample  survival  analysis  problem.  Motivation  for 
these  statistics,  results  from  small  sample  simulation  studies 
and  some  ideas  on  the  choice  of  desirable  weight  functions  from 
an  applied  point  of  view  will  be  presented  in  a  forthcoming  paper 
by  the  author.  Extensions  to  the  k-sample  problem  have  also  been 
developed  (O'Sullivan,  1986). 


APPENDIX 


Lemma  A . 


T  c 


P 


c 


unless  perhaps  if  r,  <  Tj.,, 


i 


1  or  2 ,  with  S  r,)  =0  and 


Proof. 


P 

If  Tj  =  r2  then  since  T  c  f{T  j,T  2}  and  T  j  -*  Tj  then 
P 

T  c  -»  tc.  Thus,  suppose  Tj  <  t  2.  If  tc  =  t2,  then  C  j(  t  j)  >0 

and  S  2(  rx)  =0.  Hence,  P[T  c  =  T  J  <  P[C  t(T  j)  =  0]  +  P[T  x  >  T  2] 

■»  0  by  consistency  of  C  j(T  j)  for  C  j(  Tj)  ,  and  P[T  c  =  T  2]  ■»  1.  In 

the  case  that  r2  <  t2  and  rc  =  Tj  then  either  C  2(  r 2)  >  0, 

C  t(  tj)  =0  and  S  j(  T[)  >0,  or  C  x(  Tj)  =0.  In  the  former  case, 

P[T  c  =  T  2]  <  P[S  j(T  j)  =  0]  +  P [T  2  <  T  j]  -»  0.  In  the  latter  case 

P 

if  we  also  assume  that  S  j(  t  j)  >0  then  again  S  t(T  x)  -»  S  2(  t2)  >  0 

P 

so  that  P[T  c  =  T  2]  -»  0,  and  T  e  ->  =  rc.  The  only  case 

remaining  is  that  of  r2  <  r2,  S  2(  t2)  =0  and  C  x(  t  t)  =0.  In 

general  there  is  no  guarantee  that  T  c  converges  in  this  case. 


Lemma  A. 2 


If  S  =  1  -  F  is  an  arbitrary  right  continuous  survival 
function  on  [0,  a)  and  the  variance  of  F  denoted  by  c2  is  finite, 
then 


a  2  = 


12 


!  S(u)du 
l  v 


Jo  S  ( v)  S  '(v) 


dF  ( v) 


Proof 


For  any  x  e  [0,  ®)  let 


’  x  12 

/  S  du 

<  v  i 


JO  S  ( v)  S  "(v) 


dF ( v) .  Let 


{Tn;  n=l ,  2 ,  .  .  .  }  be  a  sequence  of  real  numbers  in  [0,®)  such  that 


» 


rn<T  V  n,  {T  n}  is  increasing  and  lim  T  n  =  T  =  sup{s:  S(s)>0} 

n  ->  ® 


n  {  T  n 


;  S (u) du 


T  n  Jo  S  (V)  S  “( V) 


dF  (v) 


T  (  T 

=  I 

Jo  l  v 


S  (u)  du  dS  -1(v) 


2  T 


!  S (u) du 


S(v) 


n  ^  n 


(  /v  S(u)du) 
S  (v) 


( v )  d\ 


S(u)du)  +  2 


pr  n  v 

( ; ' 
Jo  0 


du)  S  (v)  dv 


■  t  )2  T 

!  S(u)du)  +  /  ns(v)dv2 

0  Jo 
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If  T  <  to,  then 


T  n  T  m 

lim  /  S(v)dv  2  =  T  2  S  "(T)  +  /  u  2  dF(u)  =  /u^FCu) 
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If  T  =  ®,  then 


lim  T  S(T  J  =  0  since  <r2  < 
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